Astronomy & Astrophysics manuscript no. 15087 
July 6, 2010 



©ESO2010 



O 

(N 



J. 



Herschel observations of ortho- and para-oxidaniumyl (H20'^) in 
spiral arm clouds toward Sgr B2(M) ^ 

p. Schilke,i'2 C. Comito,^ H. S. P. Muller/ E. A. Bergin,^ E. Herbst,'^ q q Lis,^ D. A. Neufeld.^o T. G. Phillips,'^ 

T. A. Bell,4 G.A. Blake,^ S. Cabrit,^^ E. Caux,'^'^ C. Ceccarelli,*^ J. Cernicharo,'^ N. R. Crockett,^ R Daniel,'^''^ 

M.-L. Dubernet/''^^ M. Emprechtinger,'^ P. Encrenaz/*' M. Gerin/'' T. F. Giesen/ J. R. Goicoechea,^ 

P R Goldsmith/^ H. Gupta/^^ C. Joblin,'^'^ D. Johnstone/^ w. D. Langeri'' W. B. Latteri^ S. D. Lord,'^ S. Maret,^ 

P G. Martin," G. J. Melnick/^ K. M. Menten,^ P Morris/^ J. A. Murphy/'^ V. Ossenkopf/2.21 l. Pagani,^^ 
C. Pearson,i3 ^ pgrault,'" R. Plume,^^ S.-L. Qin/^ M. Salez,^^ S. Schlemmer,' J. Stutzki/ N. Trappe/'^ R R S. van 
der Tak,2i C. Vastel,^-^ S. Wang,^ H. W. Yorke,'^ S. Yu,'^ ^ Erickson,^^ RW. Maiwald," J. Kooi,^ A. Karpov,'^ 

J. Zmuidzinas/ A. Boogert,"* R. Schieder,' and P. Zaal^' 

(Affiliations can be found after the references) 

Preprint online version: July 6, 2010 



in 



< 
O 

Oh 
I ' 

-(— > . 



> 

o 
i> 

p 

o 
o 



X 



ABSTRACT 

H20^ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The 
observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H20^ in most velocity intervals is 
4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the formation temperature and current 
physical temperature of the gas hosting H2O+ is discussed, but no firm conclusion is reached. In the velocity interval 0-60 km s"', an ortho/para 
ratio of below unity is found, but if this is due to an artifact of contamination by other species or real is not clear. 

Key words. ISM: abundances — ISM: molecules 



1. Introduction 

Simple di- and triatomic molecules and ions are fundamental 
constituents of interstellar chemistry which eventually leads to 
the formation of complex molecules. Many of these species 
have ground state transitions at submillimeter- and THz wave- 
lengths, and are therefore either difficult or not at all observ- 
able from the ground, yet they constitute the building blocks of 
chemistry, and are therefore fundamental to its understanding 
in various environments. Among those are spiral arm clouds, 
located in the plane of the Galactic disk, where the line-of- 
sight toward a strong continuum source passes through by 
chance. This setup allows sensitive absorption measurements, 
and the clouds have been observed by this method against the 
Sgr B2, W31c, W49, W51 and CasA millimeter continuum 
sources using molecular species s uch as CO, HCN, HCO'^, 
CS, CN, SO, a i id c-C^H? etc (e.g. [G reaves & Williams' . [1994 ; 
Tieftnink et all [19941: iGreav es & Nv man. 1996: Menten et al. . 
2010t [Gerinetali [2010r iOssenkopf et al.L [20101) . The results 
demonstrate that spiral arm clouds have low gas density and 
low excitation temperatures, and represent diffuse and translu- 
cent clouds. 

One of the best sources for these absorption studies is Sgr 
B2, located close to the Galactic center, ~ 100 pc from Sgr A* 
in projection, and one of the stronge st submillimeter sources 
in the Galaxy (e.g. [Pierce-Price et all [2000.) . The dense cores 
Sgr B2(N) and Sgr B2(M) within the cloud are at different evo- 



lutionary stages, and constitute well-studied massive star form- 
ing regions in our Galaxy. The flux ratio of the continuum be- 
tween Sgr B2(M) and Sgr B2(N) is less than unity at 1 mm and 
rises at shorte r wavelengths so that Sgr B2(M) dominates abov e 
~ 500 GHz rioldsmithetalL [l99ot [Lis & Goldsmithi [199 Ih . 
Sgr B2(M) also shows fewer mol ecular emission Hnes than 
Sgr B2(N) ( [NummeUn et al.L [1998[) . hence less confusion and 
therefore is better suited for absorption studies. The line-of-sight 
toward the Sgr B2(M) continuum will pass almost all the way to 
the center of our Galaxy, providing a more complete census in 
studying physical and chemical conditions towards the Galactic 
center clouds and all spiral arm clouds simultan eously. HIFI, the 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



Heterodyne Instrument for the Far-Infrared ( de Graau w et al. 
2010) on board the Herschel Space Observatory (Pilbratt et al. 
[2OIOI) is an ideal instrument for making these observations. 



2. Observations 

Full spectral scans of HIFI bands la, lb, and 4b towards 
SgrB2(M)(ay200() = 17''47'»20.35' and dyzooo = -28°23'03.0") 
have been carried out respectively on March 1, 2, and 5 2010, 
providing coverage of the frequency range 479 through 637 GHz 
and 1051 through 1121 GHz. 

HIFI Spectral Scans are carried out in Dual Beam Switch 
(DBS) mode, where the DBS reference beams lie approximately 
3' apart. The wide band spectrometer (WBS) is used as a back- 
end, providing a spectral resolution of 1 . 1 MHz over a 4-GHz- 
wide Intermediate Frequency (IF) band. A HIFI Spectral Scan 
consists of a number of double-sideband (DSB) spectra, tuned 
at different Local Oscillator (LO) frequencies, where the spac- 
ing between one LO setting and the next is determined by the 
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"redundancy" chosen by the observer lComito & Schilkel(l2002h . 
The molecular spectrum of Sgr B2(M) in band la and 4b has 
been scanned with a redundancy of 4, that of band lb with a 
redundancy of 8, which means that every frequency has been 
observed respectively 4 and 8 times in each sideband. Multiple 
observations of the same frequency at different LO tunings are 
necessary to separate the lower-sideband (LSB) from the upper- 
sideband (USB) spectrum. 

The data have been calibrated th rough the standa rd pipeline 
released with version 2.9 of HIPE Ott e t al.l (12010 ). and sub- 
sequently exported to CLASSQ using the HiClass task within 
HIPE. Deconvolution of the DSB data into single-sideband 
(SSB) format has been performed on CLASS. All the HIFI data 
presented here, spectral features and continuum emission, are 
deconvolved SSB spectra. Although both horizontal (H) and ver- 
tical (V) polarizations have been obtained, we will show only 
H-polarization spectra. The intensity scale is main-beam tem- 
perature, and results from applying a beam efficiency correction 
of 0.69 for band la, 0.68 for band lb, and 669 fo r band 4b 
dRoelfsemaP.. Helmich F.P.. Tevssier. P. et aUl20Toh . 



3. Spectroscopy of H2O+ 

Removal of an electron from oxidane, H2O, also known as water, 
yields oxidaniumyl, H20^. Its bond lengt hs and bond angle ar e 
slightly larger than those of H2O, see e.g. Strahan et al. ( 1986). 
Quantum-chemical calculations dWeis et al. . 1989 ) yielded a 
ground state dipole moment of ~2.4 D, considerably larger 
than in H2O. The transitions are of b-type, meaning AKa = 
AKsl mod 2. The electronic ground state changes from 'Ai in 
the neutral to ^Bi in the cation which leads to a reversal of the 
ortho and para levels with respect to water. Ka + Kc is even and 
odd for ortho- and para-li20^, respectively. The para levels do 
not show any hyperfine splitting while the ortho levels are split 
into three because of the 'H hyperfine structure. The strong lines 
have AF = A7 = AN, i.e. they do not involve a spin-flip. At 
low quantum numbers spin-flipping transitions have appreciable 
intensity. 

Further details of the spectroscopy of H20^ are discussed 
in the Appendix. Table lA. fl provides calculated rest frequencies 
for the two rotational transitions discussed in the present inves- 
tigation. Fig. lA. II show the lowest energy levels of H20^ with 
allowed transitions. 



4. Results 

Determining the opacities and thus column densities of absorp- 
tion lines is traditionally done using the line-to-continuum ra- 
tio. In the present case, this is not straightforward, because the 
orf/!o-H20^-line has hyperfine structure with closeby compo- 
nents (Fig.[T]), which distorts the simple correspondence of line- 
to-continuum ratio with column density at a given velocity, and 
also because the line background of the Sgr B2(M) core can- 
not necessarily be neglected. We therefore fitted the lines us- 
ing the XCLASSQ program, which performs an LTE fit using 
the molecular data discussed in Sect. [3j using the automated fit- 



' Continuum and Line Analysis Single-dish Software, distributed with 
the GILDAS software, see http://www.iram.fr/IRAMFR/GILDAS 

^ We made use of the myXCLASS program (https://www.astro.uni- 
koeln.de/projects/schilke/XCLASS), which accesses the CDMS 
(Mulleret al., 2001, 2005 http://www.cdms.de) and JPL iPickett et all 
ll998i[http://spec.jpl.nasa.gov I molecular data bases. 



ting routine provided by MAGI^jfl. For all velocity components, 
an excitation temperature of 2.7 K was assumed. For molecules 
that react strongly with H 2 (see the discussion in Black, 199^ 
iStauber & BrudereR l2009l) . the collisional processes in diff'use 
gas are unimportant relative to radiative excitation in controlling 
the excitation temperatures of observed transitions of species 
with high dipole moments such as H20^, since inelastic colli- 
sions with H, H2 and electrons compete with reactive collisons. 
The excitation temperature employed here may still not be en- 
tirely correct, since particularly at 1115 GHz the general FIR 
background of the Galaxy contributes to a radiation tempera- 
ture of 4.8 K even in the vicinity of the Sun, and for the spi- 
ral arm clouds one expects similar or slightly higher values 
(IWright et al.Lll991l:lPaladini et al.Ll2007h . Our analysis is not af- 
fected however if the excitation temperature is dominated by this 
radiation field, and stays significantly smaller compared with the 
upper level energies for the para- and ortho-HoO"^ ground state 
lines hvjk - 29 and 54 K, respectively, which is a well justified 
assumption. The maximum opacities of the orf/jo-H20^ line are 
about 2, so the lines are only moderately opaque. 

For para-HiO^ , only the 607 GHz line was used to perform 
the fit, since this is the strongest and least contaminated para- 
line, but it can be seen from Fig. |2]that predictions from this 
reproduce the other para lines rather well. To predict contami- 
nation, we used the fit of all species in Sgr B2(M) (Qin et al, in 
prep) as background. This is a preliminary version of the fit, and 
we cannot exclude the existence of additional contamination by 
unknown lines. Thus we estimate the error of the fit due to un- 
certainties of this nature very conservatively to be 20%, but in 
the presence of strong unknown lines it could be larger at cer- 
tain frequencies. This is particularly true for the possible D2O 
contamination of the/?flra-H20^-line at 607 GHz. The relatively 
small variation of the ortho/para ratio in the absorption cloud 
range (see below) argues against contamination by a strong un- 
known line however All para-li20^ components are optically 
thin. In the following, we make the assumption that the exci- 
tation of all upper levels can be approximated as LTE with an 
excitation temperature of 2.7 K, and that thus the ortho- and 
/:'flra-H20^ column densities can be measured by observations 
of the ground state. This assumption is reasonable for the spi- 
ral arm clouds, but most likely violated for the clouds associated 
with the Sgr B2(M) envelope (see below). 

It appears that the absorption lines of different species to- 
ward Sgr B2(M) cannot be fitted with a unique set of physical 
components of fixed velocity and velocity width. This probably 
reflects the different origins of the species in atomic, low density 
molecular and high density molecular gas with a different veloc- 
ity structure. A detailed study of the different distributions will 
have to await the complete data set of the survey. Apart from that, 
particularly in species which have hyperfine structure and are 
very abundant, that is in species which absorb at all velocities, 
the decomposition into basically Gaussian components would 
not necessarily be unique. The fit rather represents a deconvolu- 
tion of the hyperfine pattern. We therefore prefer to present the 
results as depicted in Fig.|3] as column densities/velocity interval 
and ratio as a function of velocity, as a sum over the components. 
The component of the Sgr B2(M) envelope, which is located at 
64 km s"', is most uncertain, because here the para-H20^ 607 
GHz line is most contaminated, and there the assumption of uni- 
formly low excitation temperature for all levels is most likely to 
be violated, since this is warm and dense gas with a strong FIR 
field which may dominate the excitation. 



https://www.astro.uni-koeln.de/projects/schilke/MAGIX 
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Since H20^ is expected to originate in mostly atomic gas 
at the edge of diffuse and translucent clouds, giving abundances 
relative to H2 or H does not make sense, si nce it exists ne i ther in 
pur ely atomic , nor in purely molecular gas. lMenten et alJ (1201 Ol) 
and lOin et al] (1201 Ol) give column densities of H2 of typically a 
few times 10^' cm"^ and H column densities of typically 10^° 



cm , so the average HaO^ abundance relative to the number of 
H nuclei is a few times 10"^, but could be much higher locally. 
The o/p ratio was calculated using the column densities, and has 
a mean of 4.8 + 0.5 for the spiral arm clouds, with little variation. 
We calculated the nuclear spin temperature using 



(1) 



N(para-H20^) Qp 

with AZi - 30.1 K, and Q„p the partition function of 
ort/!o-H20^//:'ara-H20^, repectively, given in the Appendix. The 
Qo^p include the rotational and nuclear spin part, and are refer- 
enced to a ground state energy of zero for both. Toward the spiral 
arm clouds, we find the mean nuclear spin temperature to be al- 
most constant at 21 ±2 K. At this temperature range, the ratio of 
the partition functions is almost equal to one, and the tempera- 
ture is mostly determined by the exponential factor 

In the velocity range from about to 60 km s"', the o/p ra- 
tio drops to unity, much below the high-temperature limit of 3. 
This does not reflect any thermal equilibrium and cannot be ex- 
plained by any known formation mechanism: at low tempera- 
tures, the o/p ratio is extremely high, since all the molecules are 
in their lowest (ortho) state, and at high temperatures the limit 
3 is reached, given by the nuclear spin statistics. We can only 
speculate about the cause of this unexpected result: it could be 
either a measurement error, or a real effect. If it is a measurement 
error, either the ortho-H20^ column density is underestimated, 
which could be caused by contamination from a strong emission 
line, or the para-H20^ column density is overestimated, which 
could be due to contamination from another absorption line, or 
the excitation temperatures deviate from the 2.7 K we assumed 
in such a way to produce this effect. The latter could e.g. be pro- 
duced by a bright FIR field at the location of the clouds. Taking 
the measured ratio at face value as the o/p ratio would imply that 
a process exists which produces ortho-H20^ and para-H20^ in 
equal amounts. 

5. Discussion 

Since there are no fast radiative transitions between para-H20^ 
and ortho-H20^, the derived nuclear spin temperature is thought 
to be determined by chemical processes, either at formation or 
afterwards, since para-ortho transformation can only occur ac- 
companied with a proton exchange reaction of one of the hydro- 
gen nuclei. The only observed kinetic temp erature estimates in 
these clouds are from' Tieftrunk et al] ( 11994 !). based on NH3, and 
show values of 35+5 K in the -100 km s~' component (which is 
believed to originate in the Galactic center), below 20 K in the 
velocity range in the spiral arm clouds, while the -10 to 20 km 
s ' component (a lso from the Galactic center) has temperatures 
exceeding 100 K ( Gardner et alifl988 ^. and is most likely shock 
heated. There is no similarity to the more or less constant nu- 
clear spin temperature of 21 K we derive for H20^ formation 
in this range, which suggests that H20^ and NH3 trace very 
different gas components: H20^ the warm outer photon domi- 
nated edge of clouds, NH3 either the shielded and cold interior, 
or hot shocked gas, which also seems to be devoid of H20^. 
This picture of H20^ formation at the edges of clouds, actu- 
ally in regions where the gas is mostly atomic, is supported by 
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Fig. 1. ortho-}i20^ , as already shown in lOssenkopf et al] ( l2010l) . 
The data are shown in black, the fit in red, in blue the hfs pattern 
is depicted, and in green the predicted contamination by other 
molecules. 








Fig. 2. para-}i20^ lines in Sgr B2, with the predicted contami- 
nation by other molecules in green, as in Fig. [1] The position of 
the D2O ground state line is indicated. 



iGerin et all (l2010l) and iNeufeld et aP (I2OIOI) based on studies 
of OH-l - and H20'^ with Herschel/UIFI. PDR models of diffuse 
clouds dLe Petit et al.. 2006) predict temperatures of 50 to 100 K 
for the formation region of H20"^ in diffuse clouds (Ay ~ 1-3) 
with densities of about 10^ and 10^ cm -^ with radiation fields of 
1-3 times ambient. 

The relationship of nuclear spin temperature and formation 
or ambient temperature needs t o be discussed in somewhat more 
detail (see iFlower et al.L l2006l for a discussion on this ratio for 
molecular hydrogen). H20^ is a very reactive ion, and partic- 
ularly it reacts exothermically with H2 to form HbO^, so the 
reaction para-li20^ -H H2 ^ ortho-li20^+ H2, which would 
equilibrate the o/p ratio to the kinetic gas temperature, might 
not be relevant. If an equivalent reaction with atomic hydrogen 
{para-H20^ -H H — > ortho-li20^+ H) could happen, is unknown. 
If it does, it most likely will equilibrate to the current gas tem- 
perature, if not, the observed o/p ratio is the one established at 
formation. How this depends on the kinetic temperature at for- 
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SgrB2(M) column densities and ratios 
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Fig. 3. Column density distribution of ortho-li20^ and 
para-H20^ (upper panel), o/p ratio (central panel) and Tnucieaispin 
distribution (lower panel). 

mation is not well understood. HiO^ is produced by the highly 
exothermic reaction of OH^ with H2, so the variables determin- 
ing the H20^ o/p ratio are 

1 . the o/p ratio of H2, 

2. how much and in which way the excess energy of the 
exothermic reaction is available for o/p conversion of H20^, 

3. the temperature of the gas at the time of formation, typically 
above 50-100 K based on PDR models. 

The latter two processes would push the o/p ratio down toward 
3:1, the high temperature value, which means toward a spin tem- 
perature higher than the 21 K we measure. The influence of the 
o/p ratio of H2 is hard to assess: if the reaction proceeds by 
way of a collision complex, then the nuclear spin of H2 will 
have an eff'ect on the nuclear spin of product H20^, but not if 
the reaction proceeds through direct atom trans fer. The mea- 
sured o/p of H2 in diffuse clouds is about unity (Savag e et al.L 
[1977; Rachfordet al., 2009). Clearly, this is an interesting area 
of molecular physics that needs further study. From our obser- 
vations, it seems that we see an excess of ortho-H20^ relative 
to what one would expect based on the formation temperature 
and available reaction energy in the spiral arm clouds and, if 
the measurement of an ortho-H20^ Ipara-HiO^ ratio in the to 
60 km s ' region is real, an excess of para-H20^ there. This 
velocity range nominally is assigned to the Sagittarius/Scutum 
arms ([^llee, 2008), outside of the Galactic center, where no ex- 
otic conditions are expected. However, this velocity range is also 
bracketed by gas local to Sgr B2 at to 10 km s"' and around 60 
km s"', so it could represent diffuse gas belonging to this com- 
plex, which could be the cause of unusual excitation or chemical 
conditions, e.g. due to shocks. In the lower velo city range over- 
lapping with this (-10 to 20 km s~'). lLis etaP (12010 ) also find 
wa ter with an o/p ra tio of 3, indicating high temperatures. 

JLis et al.l (l2010l) find an average spin temperature of 27 K 
for water toward the the spiral arms lines of sight. H2O in these 
clouds is formed in the gas phase, through dissociative recombi- 
nation of H30^, in a region where the gas is at about this tem- 
perature. The correspondence between physical temperature and 
spin temperature may be more easily traced by the more stable 
water molecule, although in general the contribution of grain sur- 
face chemistry for water complicates the issue (see discussion in 
iLis et alj|2010h . From this study it is clear that by determining 
the ortho/para-raiio of H20^ (and, by proxy, from other simple 



hydrides) one can learn a lot about the formation processes, but 
also that many fundamental physical and chemical processes are 
still not fully understood. We can look forward to the wealth of 
data HIFI will bring ! 
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Table A.l. Quantum numbers of rotational transitions of HaO^ 
described in the present work, calculated frequencies (MHz) 
with uncertainties in parentheses"; lower state energies £10 (K) 
and Einstein A-values ( 1 0"^ s" 1 ) 



para 



Elk (K) 



N' , , 
K'K. 



■N"„ „ 

K'K" 



J' - J" 



F' - F" 



frequency 



li 



Im, para 



1.5- 
1.5- 
0.5- 
0.5- 



0.5 
1.5 
0.5 
1.5 



604678.6 (25) 
607227.3 (19) 
631724.1(37) 
634272.9 (24) 



0.005 
0.000 
0.005 
0.000 



111 - Ooo, ortho 



1.3 

6.2 
5.6 
2.8 



1.5- 


-0.5 


1.5- 


-0.5 


1115150.75(85) 


0.122 


17.1 


1.5- 


-0.5 


0.5- 


-0.5 


1115186.18(81) 


0.122 


27.5 


1.5- 


-0.5 


2.5- 


-1.5 


1115204.15(82) 


0.000 


31.0 


1.5- 


-0.5 


1.5- 


- 1.5 


1115262.90(82) 


0.000 


13.9 


1.5- 


-0.5 


0.5- 


-1.5 


1115298.33(87) 


0.000 


3.5 


0.5- 


-0.5 


0.5- 


-0.5 


1139541.54(103) 


0.122 


3.7 


0.5- 


-0.5 


1.5- 


-0.5 


1139560.58(94) 


0.122 


14.8 


0.5- 


-0.5 


0.5- 


- 1.5 


1139653.69(94) 


0.000 


29.4 


0.5- 


-0.5 


1.5- 


- 1.5 


1139672.73(103) 


0.000 


18.3 



" Numbers in parentheses are 1 a uncertainties in units of the least 
significant figures. These values should be viewed with some caution, 
in particular for the para transition, see text. 

' F is redundant for para transitions; F = J may be assumed. The 
lowest para state is 30.01 K above the lowest ortho state. 




orttio 




Fig. A.l. Detail of the energy level diagram of H20^. Hyperfine 
splitting has been omitted for the orf/10-levels. Rotational level 
assignments Nk„k, are given below the levels, fine structure 
level assignments J to the side. Magenta arrows mark transi- 
tions observed in the course of the present investigation. All 
other transitions shown can be observed with HIFI. The thick- 
ness of the arrows indicates the relative strengths of the transi- 
tions and the numbers the approximate frequencies. Frequencies 
of the weaker components are given in parentheses. The only 
transitions connecting to levels with N - I which are not shown 
are 221 - 1 10 and 220 -In near 2.85 and 3.0 THz which are not 
observable with HIFI but with PACS. 



Appendix A: The spectroscopy of H2O+ 

The rotational spectrum of o xidaniumyl was measured by laser 
magnetic resonance (LMR) ("Straha n et al.L Il986t iMtirtz et all 
[1998); further infrared a nd electronic spect ral measurements 
have been summarized in IZheng et all (l2008h . Observations of 
the Nk„Kc = 111 - Ooo, J - 1.5 - 0. 5 fine structure compon ent 
near 1115 GHz with Herschel/mPl (lOssenkopf et all 1201 Oh as 
well as subsequent observations raised the issue which of the two 
sets of spectroscopic parameters from LMR measurements pro- 
vide more reliable frequency predictions. Latest observations as 
well as reint erpretations of older ones favor the parameters from 
iMiirtz et aP ([1998) even though there seem to be small discrep- 
ancies of order of +5 MHz or 1.35 km s' between various obser- 
vations for the specific transition mentioned above. Observations 
carried out toward Sgr B2(M) or Orion KL for the HEXOS pro- 
gram are probably less suited to derive rest frequencies than cer- 
tain other observations. Therefore, a preli minary catalog entry 
has been constructed for the CDMS catalog (Miiller et al., 2001, 
[2005); the final entry is intended to be a common CDMS and 
JPL JPickett et all 11998]) catalog entry. 

All infrared data and all ground state combination differ- 
en ces (GSCDs ) deriv ed from electronic spectra as summarized 
in lZhen^ et al.l (i2008h were used in the fit as long as they were 
deemed reliable. These data are uncertain to between 0.005 and 
0.030 cm-' or between 150 and 900 MHz. iMiirtz et alJ (Il998l) 
provided for their measured data extrapolated zero-field frequen- 
cies as well as residuals between observed and calculated fre- 
quencies along with uncertainties. From these data weighted av- 
erages of the hypothetical experimental zero-field frequencies 
and of their uncertainties were derived; these uncer tainties were 
of o rder of 2 MHz with a considerable scatter S trahan et alJ 
(11986) do n ot give sufficient data for this purpose. However, 
IMiirtz etal] (11998) published calculated frequencies for the 1 1 1 - 
Ooo transition. Because of the importance of this transition for 



the astronomical observation as well as for the fit, these cal- 
culated frequencies were also used in the fit with presumable 
uncertainties corresponding to 2 MHz. The quantum numbers, 
calculated frequencies, and uncertainties for the two observed 
rotational transitions are given in Table lA.ll Even though the 
LMR data fit well within their uncertainties, the calculated fre- 
quencies should be viewed with some caution because it is not 
clear how reliable the zero-field extrapolation is. Moreover, the 
large centrifugal distortion effects affecting the spectra of this 
ion require additional caution with respect to any extrapolation. 
It is worthwhile mentioning that the In- Ooo transition frequen- 
cies in that table are essentially identical to the ones given in 
Mu rtz et al.1 ([1998). The situation is diff'erent for the lu) - loi 
transition. The J - 1.5-1.5 fine structure splitting derived from 
term values given in Table V of Miirtz et al. (1998) differs from 
the value in Table [ATI by less than 4 MHz whereas the remain- 
ing fine structure intervals differ by up to 80 MHz. On t he other 
hand, calculations directly from the IMiirtz et al] (1 19981) param- 
eters differ from values in Table [ATI by less than 10 MHz. As 
can be seen in Fig. lA.ll the ground state para level loi of H20^ 
is -30.13 and 30.01 K above the Ooo level for the J = 0.5 and 
1.5 fine structure components respectively. Both para and ortho 
ground state levels are split into 2 because of the fine and hyper 
fine structure splitting, respectively. The quantum numbers are 
0.5 and 1.5 in both cases, giving g = 2 and 4, respectively, and 
hence Q x 6 at low temperatures if ortho and para states are 
treated independently. This 1 : 1 ratio for Q at low temperatures 
approaches 3 : 1 at room temperarture; it is about 1 .5 : 1 and 
2 : 1 at ~40 and 75 K, respectively. Table lA.2] gives selected par- 
tition function values for para- and ortho-H20^, assuming they 
are completely non-interacting species. 

The mixing of ortho and para states can be mediated by 
terms such the off-diagonal electron spin-hydrogen nuclear spin 
coupling term Tah or the offOdiagonal hydrogen nuclear spin 
coupling term Cat + Cta- The radical NH2 and PH2 are isoelec- 
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Table A.2. 



T 


Qipara) 


Qiortho) 


300.0 


109.1787 


296.3567 


225.0 


73.4887 


192.9372 


150.0 


43.0579 


105.7512 


75.0 


19.0770 


38.3647 


37.5 


10.4751 


14.5910 


18.750 


7.3879 


7.3124 


9.375 


6.2329 


6.0584 


5.000 


5.9680 


5.9982 


2.725 


5.9123 


5.9961 



tronic and isovalent to H20^. Model calculations have shown the 
largest perturbations to occur between the loi and In levels, but 
they are with less than 5 (Gendriesch et al., 2001) and less than 
3 kHz (IMargules et al.i, 2002), respectively, rather small, maybe 
even negligible. Model calculations suggest that perturbations of 
the two rotational transitions described in the present study are 
less than 3 kHz. Slightly larger perturbations may occur at higher 
quantum numbers. 



